Exponential distributions of collective flow-event properties in viscous liquid dynamics 
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We study the statistics of flow events in the inherent dynamics in supercooled two- and three-dimensional 
binary Lennard-Jones liquids. Distributions of changes of the collective quantities energy, pressure and shear 
stress become exponential at low temperatures, as does that of the event "size" S = ^ di 2 . We show how the 
S-distribution controls the others, while itself following from exponential tails in the distributions of (1) single 
particle displacements d, involving a Lindemann-like length g(l and (2) the number of active particles (with 
d > d L ). 



Many complex systems, including turbulent flows [1], plas- 
tically deforming crystals [2], and financial markets |2|], are 
characterized by intermittent, stochastic dynamics. It is im- 
portant to characterize the statistics of the individual dynam- 
ical events, called "increments" in formal stochastic models, 
"returns" in finance, and "avalanches" in driven systems [4]. 
In particular non-Gaussian effects can have crucial signifi- 
cance, as extreme events can dominate the dynamics. For 
example, in finance the fact that distributions of returns are 
known to have "fat" (i.e. power-law) tails — in opposition 
to assumptions made by standard models J5[] — has been pro- 
posed as part of the reason for the failure to understand mar- 
ket crashes JfJ. In Gaussian statistics, values more than five 
standard deviations from the mean account for less than one 
millionth of the distribution. For a symmetric exponential 
distribution, on the other hand, they account for about one 
thousandth. For the fat-tailed case of a Cauchy distribution 
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the equivalent fraction (in terms of the stan- 



dard deviation of the best-fit Gaussian) is nearly one tenth. 
Examples of exponential distributions (which may be consid- 
ered intermediate between the Gaussian and power-law cases) 
include bursts of protein production from gene expression [7] 
and even returns in financial time series (at least for not-too 
large values) 0]. 

Rearrangements of molecules known as "flow events" are 
the fundamental processes giving rise to structural relaxation 
and flow in disordered systems ranging from the equilibrium 
case of highly viscous liquids jSJ iKJ QJJ] to non-equilibrium 
systems such as glasses J3l, granular materials [12] and foams 
II 1 311 - For viscous (or "glass-forming") liquids, it has been ac- 
cepted since Goldstein's 1969 paper [9] that the dynamics may 
be understood in terms of a division into fast vibrational mo- 
tion around a particular energy minimum, and relatively rare 
transitions between neighboring minima. It is the latter that 
are responsible for the slow dynamics ifloL 14], thus a de- 
tailed understanding of their nature is essential, particularly 
since many theoretical approaches start by making assump- 
tions about the flow events lfl5i[l6l[l7ll . Computer simulations 
provide the means to go beyond assumptions, and in recent 
years have provided much insight into the kinds of molecu- 
lar motions that occur in a supercooled liquid. Most workers 
have concentrated on particle displacements. But of more di- 
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FIG. 1: (Color online) Distributions of A N E, A N p, A N a B for T = 
0.33 (2D); dashed lines show exponential fits (without transforming 
by taking logarithms) to values above 0.3, with decay lengths 0.160, 
0.162, 0.284 for E, p, a s , respectively. Inset, distribution of Aa s 
at T=1.5 and 0.33, now normalized to have unit standard deviation, 
and compared with a standard normal distribution. Non-Gaussian 
parameters (NGP) are indicated for each distribution. 



rect relevance to experiments are collective properties of the 
dynamics — changes in quantities such as potential energy E, 
pressure p, and shear stress a s . Knowledge of their distri- 
butions, and how those depend on temperature, is sure to be 
relevant for understanding the slowing down of dynamics in 
viscous liquids. Vogel et al. IU8I1 studied the relation between 
particle displacements and the size of energy changes dur- 
ing flow events, and while they reported exponential tails in 
both, they did not examine their relationship in detail, while 
Schr0der et al. also reported such tails in flow event displace- 
ments [ 10]. Exponential distributions of forces in a simulated 
Lennard-Jones fluid were observed already in 1987 11911 . al- 
though their relevance to flow event properties is unclear. It is 
now becoming possible to study individual particle displace- 
ments experimentally, as Schall et al have done for a colloidal 
system [20]; as yet the analysis is limited to a few events. The 
exponential tails recently reported for the van Hove correla- 
tion function 112 ill reflect the cumulative effect of many flow 
events, whereas this work focuses on statistics on single flow 
events. 

In this Letter we present simulation data from two- and 
three-dimensional (2D and 3D) binary Lennard-Jones (BLJ) 
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model fluids, brought to a viscous state by cooling. By study- 
ing the so-called "inherent dynamics" II 1 C XI — the trajectory ob- 
tained by mapping configurations to local energy minima — 
we identify flow events and study their statistics. Our main re- 
sults are (1) at low temperatures T the distributions of changes 
of E, p and a s are exponential, in contrast to high tem- 
peratures where they are basically Gaussian; (2) the sum of 
squared displacements S, a geometrical measure of the size of 
an event, is the controlling quantity in the sense that it also has 
an exponential distribution (ED) at low temperatures, but for 
fixed S values E, p and <j s have Gaussian distributions. The 
mean event size decreases with decreasing temperature; (3) 
the ED of S can be traced to the existence of an exponential 
tail in the distribution of particle displacements during events, 
characterized by a Lindemann-like length scale, which defines 
"active" particles. The number of active particles is broadly 
distributed at high T, typically comprising a large fraction of 
the particles in the system, whereas at low T it is also ex- 
ponentially distributed, with a mean of a few tens of particles. 
This crossover coincides with the increasing relevance of min- 
ima transitions to the dynamics at lower temperature. Indeed, 
at high temperature a transition occurs almost every time step 
and their relevance to dynamics is clearly minimal, whereas 
at lower temperatures events occur one at a time in a localized 
part of the sample. We leave analysis of the waiting times and 
correlations between events to later work. 

The parameters for the BLJ potential are (where e and a are 
the energy and length scales for interactions between large (L- 
) and small (S-) particles) ell = 1, £ss = 0.5, cls = 1-5, 
<7ll = 1, agg = 0.88, gls — 0.8. All particles have the 
same mass m — 1. These parameters are identical to those of 
the BLJ introduced by Kob and Andersen 1E2I1 . The potential 
was truncated using an interpolating polynomial between 2.4 
a a /3 and 2.7 cr a p (a,/3 s {L, S}). All results reported here 
are from constant volume simulations with periodic bound- 
ary conditions, with N p =700 (1372) particles in 2D (3D), of 
which 60% (80%) were of type L, at a density of 1.2cr^ 
(1.2ct^"!), using a time step of 0.01 <JLL\Jm/ ^ll (0.005 at 
T >1.0). The system size was chosen to allow a large range 
of event sizes. From now on, all quantities will be reported 
in the "natural" units defined by cll, &ll and m. Two-step 
relaxation, a signature of landscape-influenced dynamics, first 
appears at T ~ 0.50 (2D) and T = 0.60 (3D). 

To identify flow events we carried out Stillinger's proce- 
dure llzjl of quenching configurations to the "nearest" local 
minimum. For T > 0.5, we quenched every time step. For 
lower T we quenched every tenth time step, and if a change 
of minimum was observed, the simulation was "backtracked" 
and quenched at each of the intervening time steps. We de- 
tected events using the sum of squared particle displacements 
S = J^i dj, where di is the magnitude of the displacement 
of the zth particle between successive inherent structures: A 
value greater than 10~ 3 is sufficient to distinguish a genuine 
change of minimum from numerical noise. This criterion is 
consistent with one based on changes in inherent energy or 
stress. Care was exercised in the minimization process II 3 311 . 
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FIG. 2: (Color online) (a) Distributions of S for events for different 
temperatures, and fits to the tails (Gaussian for T=1.0, 1.5, other- 
wise exponential). Structural relaxation times r a determined from 
the self-intermediate scattering function are also indicated, (b) Left 
panel, normalized change in shear stress for event size S in different 
narrow intervals at T — 0.33. Right panels, variance of (absolute) 
changes of (top) E and (bottom) <r s as a function of S. For the range 
of S relevant for lower T, the variances are almost linear in S. This 
linear relation is independent of T for cr„, while the slope increases 
with T for E. 



In the following we place somewhat more emphasis on the 2D 
data because better statistics were obtained; due to a larger 
number of particles in 3D (though the linear size is more than 
a factor of two smaller) and the larger number of neighbors per 
particle, the minimization is more time consuming (by about a 
factor of six). Thus 10-20 x 10 3 events per temperature were 
obtained in 3D compared to 10 5 in 2D. 

Figure Q] shows P(A N X) as a function of |A W JT|, where 
X is E, p or a s and the A N denotes the change normal- 
ized by the r.m.s. equilibrium fluctuations {-J {{X — (X)) 2 )) 
calculated from the time series of inherent states (T=0.33). 
The tails are very close to exponential for |A W X| > 0.25. 
Though the distributions differ at small values, the tails for p 
and E are very similar, with the same decay length A ^0.16 
(possibly a reflection of the strong-pressure correlations re- 
cently reported for Van der Waals systems [24]), while that 
for shear is larger, A ~ 0.28. To quantify how close to 
Gaussian a distribution is we use a non-Gaussian parameter 
NGP= ((fa) 4 )/3((((5x) 2 }) 2 - 1 for any quantity x. This is 
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FIG. 3: (Color online) Left panel, distributions of L-particle dis- 
placements during events at four temperatures, and exponential fits 
of tails (d >0.3), with A denoting decay length. S-particle displace- 
ments are shown for the lowest temperature (where the difference 
is greatest). Right panel, distribution of Na, the number of active 
particles, for the same temperatures, with fits for the three lowest 
temperatures. A visualization of the 35 active particles in an event at 
T=0.33 is shown in the upper part. 



zero for a Gaussian distribution and unity for a (symmetric) 
ED. The inset of the figure compares distributions of shear 
stress changes at the highest and lowest temperatures, now 
normalized to unit width. That from T=1.5 is clearly more 
Gaussian, albeit with an exponential tail a few standard devi- 
ations from the mean. 

By symmetry, changes of shear stress are uncorrelated with 
those of energy and pressure, meaning that none of these 
quantities can be considered a meaningful measure of the 
"size" of an event. A natural measure of the size of an event 
is the quantity used to detect them, S = J2i d%- The S dis- 
tributions are shown in Fig. Ufa). Going from high to low 
temperature, there is a dramatic reduction in the average size, 
while the shape of the distribution becomes increasingly expo- 
nential |34ll . Examining the distributions of AX for a narrow 
range of S, we find they are Gaussian even at low T when 
the full distributions are exponential (Fig. |2b)). For not too 
large values of S, {S < 10), the variance for these fixed-5 
distributions rises approximately linearly with S (right panels 
of Fig. [2b)). This shows that for a given S, AX is essentially 
a sum of random contributions, whose number is proportional 
to S, the size of the event. This relationship is independent of 
T for Aa s (lower right panel), but this holds only for small S 
for Ap and AE (upper right panel). It is not clear why this is. 
In all cases the variance tends to saturate as S exceeds 20 (50 
for^atT = 1.0). 

The above can be made mathematically explicit by writing 
P{AX) = J °° P(AX\S)P(S)dS. At low T if P(S) = 
(1/(5)) exp(— S/(S)) and the conditional probability 
P(AX\S) is a Gaussian with variance ax S, then integration 
gives P(AX) = l/(2a x (S}) exp(-2|AX|/ \/2a x (S)). 
Linear fits to the S < 10 data for T=0.33 gives for ax 



the values 0.32, 1.5 xl0~ 3 and 5.8 xlO" 3 for E, p and 
<7 S , respectively. The decay lengths of the corresponding 
EDs should be y / a x (S)/2. For T=0.33, (5)= 1.83, we 
get values 0.54, 0.012 and 0.023. After normalizing by the 
r.m.s. fluctuations as was done in Fig.Q]these become 0.16, 
0.15, 0.26, in reasonable agreement with the decay lengths 
determined in Fig. Q] 

To understand the distributions of S we consider the distri- 
butions of individual (large) particle displacements d in flow 
events, shown in Fig. [3] As found in Refs. fiol and ITii clear 
exponential tails appear for d larger than about 0.2. The de- 
cay lengths vary surprisingly little over the simulated temper- 
ature range, 0.08-0.15. Most of the variation is in the number 
of particles in the tail region (the average number per event). 
We find it useful to take the length scale ~ 0.1 suggested by 
the decay length seriously, and define the "active" particles as 
those with d > di=0. 1 (independent of T for simplicity). 

Distributions of Na, the number of active particles, are 
shown in the right panel of Fig. [3] There is a striking shift 
in both the shape and the mean value from high to low T: 
A typical event at high T involves most of the system; in- 
deed, larger events are more probable than smaller ones. At 
low T the distribution becomes exponential with a mean of 
order a few tens of particles. In this regime, the events are 
relatively localized though not necessarily compact (see the 
example in Fig.^-string-like spatial correlation is present in 
varying degrees in small and intermediate-sized events, while 
the largest events tend to have a more compact structure. The 
nature of spatial correlations/structure of flow events has been 
investigated by different authors. While Schr0der et al. found 
string-like correlations for transitions between minima, Ap- 
pignanesi et al. used the distance matrix method to identify 
transitions between so-called meta-basins 12511 . Their analy- 
sis highlights events involving a large fraction of the system, 
which was relatively small (150 particles). In fact their data 
(Fig. 2 of Ref. 1251) is consistent with an exponential tail with 
characteristic size of order 10-20 particles; the "democratic 
events" would simply be the largest events in the tail. On the 
other hand we observe string-like features in smaller events, 
but not so much in larger ones. Thus our results encompass 
both those of Schr0der et al. and Appignanesi et al. 

Fig. @] shows data for a 3D system. The basic features are 
similar to the 2D case, and in particular the S distribution con- 
trols the others. Also shown in the third panel Fig. [4] is the 
distribution of Na for a smaller system, iVp=110, which is 
clearly cut off by the system size — this could well explain a 
factor of four difference in relaxation time between the two 
sizes that we observe. 

The notation d^ is meant to suggest a Lindemann-like inter- 
pretation. Several authors have emphasized the importance of 
such a length scale in the context of mechanical instabilities 
associated with structural relaxation in liquids Jl7, 2(| 271. 
also in experiments H20fl . Its appearance in the present con- 
text can be interpreted as that events are in a loose sense "dis- 
cretized" in units of c?l: each event involves some number of 
particles each of which is displaced some number of di, units. 



4 



I 



. T=0.60, T a =28 
. T=0.48, T o ,=520 
» T=0.42, x =25000 



...«T=0.48, N=110 



X 

*♦ "■ < • «♦••» 

*♦ ". 5. . ■*. 



•i* •••• 



10 20 30 40 0.2 0.4 0.6 0.8 1 50 100 150 

S d N A 



FIG. 4: (Color online)Distributions of S, d (L-particles) and Na for 
a 3D system at temperatures T=0.60, 0.48 and 0.42. P(Na) for a 
system with Np=l 10 at T = 0.48 is also shown. 



Exponential distributions may be interpreted as the statement 
that there is no preferred number of basic units (this is analo- 
gous to a simple model for polymer size distributions, where 
the probability to attach a new monomer does not depend on 
the current length of the chain ll28lo . 

The cross-over to exponential distributions is associated 
with one to well-defined localized events, whose mean size 
decreases with T. The decrease indicates increasing material 
stability. A similar result was found by Fabricius and Stari- 
olo when perturbing equilibrium configurations and compar- 
ing the corresponding inherent states 12911 . If event sizes de- 
crease with decreasing temperature, this presumably has con- 
sequences for relaxation times, because a given change takes 
more events, although decreasing size (in particular AE) also 
suggests that energy barriers for individual events also de- 
crease (we will study energy barriers in upcoming work). On 
the other hand, when there is a distribution of barriers, there 
should be increased tendency for repeated reverse-transitions. 
This could be examined by studying correlations between 
events, which we have ignored here, but will address in up- 
coming work. One way to account for forward-backward 
correlations, proposed by Doliwa and Heuer l30ll . involves 
grouping minima into larger structures called metabasins and 
study transitions between these. Their analysis indicates 
growing effective barriers as T decreases, consistent with 
non-Arrhenius slowing down of dynamics. Our analysis fo- 
cuses on the more fundamental units of dynamics, individ- 
ual minima-transitions, because these are (in principle) unam- 
biguously defined and it is worth understanding their statis- 
tics in detail before attempting to coarse-grain. The analysis 
suggests the somewhat paradoxical result that individual event 
barriers decrease as T does. This emphasizes even more the 
role of correlations. 

Center for viscous liquid dynamics "Glass And Time" is 
sponsored by The Danish National Research Foundation. 
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